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Abstract. - The size segregation of granular materials in a vibrating container is investigated 
using Molecular Dynamics. We find that the rising of larger particles is accompanied by the 
existence of convection cells even in the case of the lowest possible frequencies. The convection 
can, however, also be triggered by the larger particle itself. The possibility of rising through 
this mechanism strongly depends on the depth of the larger particle. 



One of the most puzzling phenomena encountered in granular matter is size segregation: 
When a mixture of grains of the same material (equal density) but different size is shaken 
in a container the larger particles rise to the top. This effect has been extensively studied 
experimentally 0, ^ and has much importance in numerous industrial processes Q . Recently 
this so called "Brazil nut effect" has attracted much interest among physicists 

Size segregation inevitably seems to contradict equilibrium statistical mechanics since the 
density of the overall packing increases with polydispersivity and so gravity makes situations 
with larger particles on the bottom energetically more favourable. Rosato et al.[^ proposed 
a Monte Carlo algorithm and put forward a kinetic argument to explain segregation using 
the fact that smaller particles are more mobile. In the same year Haff and WernerQ did 
Molecular Dynamics simulations of rather small systems and claimed that segregation was 
essentially a consequence of solid friction and the rotation of the particles. JuUien et al.[|| 
used a piling technique which is non-stochastic as compared to the one of ref. 6 and found a 
critical ratio TZ for the radia of spherical particles below which no segregation occurs. Based on 
these ideas Duran et al.[^ formulated a geometrical theory for segregation in which the small 
particles glide down along the surface of the larger particle, the critical ratio TZ of radia being 
between continuous and discontinuous gliding. They also claimed experimental evidence for 
two types of dynamics and visualized the discontinuous ascent of the larger particle through 
stroboscopic photos. JuUicn et al.[^ reproduced the discontinuous dynamics by including 
horizontal random fluctuations into their model. 

Parallel to these local theories there has been the "convection connection": It is known 
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experimentally[|rT| and numerically |l3| that shaken assemblies of spheres form convection 
rolls which are attached to the walls of the container. For weak shaking the convection rolls 
only appear on the surface. Knight et al.[0 showed experimentally that segregation was due 
to this convection and the fact that larger particles have a harder time entering again into the 
bulk once they are on the surface. They also verified an exponential decay of the convection 
strength as function of depth for weak shaking Jl5|. Duran et al.[^ verified segregation due 
to convection in two dimensions for strong shaking and claimed that above mentioned local 
mechanisms are at work at weak shaking. Their stroboscopic pictures showed convection cells 
above the particles. 

Here we present large scale Molecular Dynamics simulations showing that also for weak 
shaking convection is responsible for segregation but in a more intricate way: Under certain 
conditions the larger particle is able to pull down the convection rolls due to the more efficient 
momentum transfer and then rises within its flow. Because of the exponential decay of 
convection with depth the ability to rise critically depends on the vertical position of the 
larger particles. 

We used the classical Gear predictor-corrector Molecular Dynamics [O] algorithm and the 



contact forces proposed by Cundall and Strack|18 : Two grains i and j at positions f,; and 



r 



interact if the distance between the center points is smaller than the sum of their radii 
I — I < Ri + Rj. For this case the force between particles i and j moving with the 
velocities fi and fj and rotating with angular velocities fli and flj is given by 

N "^i " S f ^ ^ ^3 

with the normal and shear forces being 

FiJ ^ Y- {R, + i?^. - |r, - r, I ) + 7^ ■ m^/^ • |rl - | (2) 



= sign (|t;;7'|) min {75 • ■ | , • l^^^l} . (3) 



The relative velocity v[.^^ between i and j and the effective mass ml^A^ of the particles i and 



j are defined as 

eff 



{h-r'j) + R,-n, + R, -dj (4) 



rrii ■ m 



The resulting momenta Mi and Mj acting upon the particles are 

= • i?. , M, = -i^f • R, . (6) 

Eq. (0) takes into account that the particles slide on each other for the case that the inequality 
I F[j 1 < I Ffj I holds, otherwise they feel some viscous friction. 

Throughout our simulation we used the parameters F = 3 • 10^ g/sec^ (Young modulus), 
7iv = 100 sec~^, 7s = 1 sec~^ (phenomenological normal and shear friction coefficients) 
and n = 0.5 (Coulomb friction parameter). We considered N = 950 particles with radii 
uniformly distributed in the interval R^ G [0.85, 1.15] cm and with masses m, = 2'kR^p, 
p = 1 g/cm^. The particles are put into a two-dimensional box having walls made of particles 
with the same material characteristics as the grains and which vibrates vertically according 
to ybox = A ■ sin(27r/t) with A — 2 cm. Gravity acts in negative j/-direction Ff^ = —mi g, 
g = 981cm/sec^. The time step for the numerical integration of the Gear predictor-corrector 
scheme of 5th order was At = 5 ■ 10~^sec. 
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We investigate segregation and convection behaviour as function of the vibration frequency 
/ in two different systems, either all particles are small {Ri € [0.85, 1.15] cm) or we add one 
single big particle of i?i = 4 cm located in the centre of the box close to the bottom. In order 
to investigate closer what happens at the onset of segregation we keep all the other parameters 
fixed. Fig. ^ shows the convection cells without the larger particle (left) and with the larger 
particle (right) for three different frequencies / = 2.6 sec~^ (upper figure), / = 2.8 sec~^ 
(central figure) and f — 3 sec^^ (lower figure). For f — 2.8 sec~^ the convection cells of both 
systems, with and without the big particle, differ significantly while they are quite similar in 
the other cases. This indicates that there is a certain frequency interval where the presence of 
the big sphere triggers the onset of convection which finally leads to segregation. Indeed we 
find that convection is always present when segregation happens. It is important to note that 
if one is close enough to the onset of segregation just putting the larger particles one row lower 
can entirely supress the effect of segregation. This dependence of segregation on the height is 
quite strong and has not been discussed in the literature. By changing the frequency / very 
slowly and measuring the convection fiow through a plane at a certain height we observed that 
the transition from the fluctuation regime to the convection regime is very sharp and differs 
between different runs for less than A/ — 0.05 sec~^. Moreover when increasing the frequency 
the transition occurs almost exactly at the same frequency as when decreasing the frequency, 
i.e. there is no hysteresis. This is in contrast to the case of viscous fluids which tuned out 
of equilibrium by applying a temperature gradient have nonequilibrium dissipative structures, 
e.g. Benard convection, with a characteristic hysteresis. 

The triggering of convection cells by the big particle is investigated more quantitatively by 
calculating the convective flux $ defined as the sum of material (mass) fiow in the centre of the 
box jtop and the flow close to the walls jtot by considering that these flows have opposite signs 
as illustrated in Fig. ^j. The flows jtop and jbot are defined by adding the number of particles 
which move in one direction minus the ones moving in the opposite direction and is measured 
in particles per cycle time of the box. In fact we measure for each particle if the positions at 
subsequent nodes of the vibration are on different sides of a height line, where the height of 
the box was divided into 80 height lines between the bottom of the box and the surface of the 
packing. Fig. || shows the convective flux $ through planes at different heights D for both 
systems and for the three different frequencies. For / = 2.6 sec~^ we find almost no directed 
flow but only fluctuations. For / = 3 sec~^ both systems behave similarly, as could also be 
observed in Fig. ^ For / = 2.8 sec~^, however, the convection cells clearly extend deeper 
due to the existence of the larger particle. Apparently the larger particle is able to pull the 
convection cells down. 

This effect can be explained by the fact that in the region around the large particle the 
accelerations are higher since the momentum is transferred with less dissipative loss through 
the larger particle than through a corresponding pile of smaller particles of the same volume. 
We measured the sum of the absolute values of the forces of all the particles in a region around 
the position of the large particle and averaged it over time. The region was ring shaped 
with the inner border of radius 4.5 cm and the outer border of radius 8 cm. At the onset of 
segregation (/ = 2.8 sec~^) the average force of the small particles around the large particle 
is about 15% larger than that of the small particles in the same region in a system containing 
no large particle. This effect is strongest in the lower part of the ring shaped region. For 
/ = 2.6 sec~^ the difference is only 5%. Therefore the accelerations in the region around a 
large particle are larger than if no particle would be present. We believe that this increase in 
high frequency oscillations is responsible for pulling the convection rolls down. It is, however, 
interesting to note that the granular temperatures in this region is roughly the same in the 
two systems. 
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Fig. 1. - Convection rolls in systems without (left) and with the larger particle (right) for three 
different shaking frequencies, / = 2.6 sec~^ (upper figure), / — 2.8 sec~^ (central figure) and 
f — 3 sec~^ (lower figure). In the central case the big particle triggers convection rolls. The pictures 
were obtained by averaging over 50 shaking periods of the box. 



One can see from Fig. ^ that the convection cells decay very sharply in strength but that 
even in the deep regions some essentially horizontal motion occurs. This is reminiscent of the 
stroboscopic pictures of [ p^ implying that even in the low acceleration regime some particles 
move inward horizontally. Within our framework this motion could be interpreted as the 
exponentially weak tail of the convection rolls. 
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Fig. 2. - Schematic diagramm illustrating how the material flow jbot close by the walls and jtop in the 
center were added to obtain the convective flow <1>. 



Next we investigated the dependence of the onset of convection on the ratio of radia TZ. 
Note that in our case TZ = Ri/lcm because the mean radius of the small particles is 1cm. For 
/ = 3.2 sec"^ we observed a big particle of radius TZ = 4.0 moving up immediately. We also 
studied the cases TZ — 3.5, TZ — 3.0, TZ — 2.5 and TZ ~ 2.0. In these cases the large particle 
remains a certain waiting time on the bottom before it suddenly moves up quite rapidly. Fig. 
^ shows a typical evolution of the vertical position of the big particle with radius TZ — 2.0 
and TZ = 3.0. The waiting times do not noticeably depend on TZ being of the order of 30 sec. 
Once the large particle comes to the top it performs an oscillating motion going up and down 
(whale effect) that has also been observed experimentally [p9[ . This motion seems due to the 
convection rolls: In the case where TZ is smaller the oscillating motion is more regular because 
the larger particle has less difficulty reentering into the bulk from the surface and to follow the 
convective motion. For larger TZ the larger particle has more difficulty reentering thus leading 
to a more erratic horizontal motion as seen in Fig. ^. Particles with smaller TZ also seem to 
dip deeper into the bulk showing that the convection cells can move them more efficiently. 

We have shown with large numerical calculations of granular media in a vibrating box that 
in two dimensions segregation is intimately connected to convection. The larger particles, sur- 
rounded by a region of higher acceleration, loosen the material thus deepening the penetration 
of the convection rolls from the top into the medium. After some waiting time the larger 
particles are caught by the lower part of the rolls and pulled up. This triggering effect is only 
relevant close to the characteristic onset frequency of segregation which is sharply defined and 
strongly dependent on the initial depth of the large particle. Once the large particle is on the 
top it periodically goes up and down (whale effect) driven by the convection cells. 

In addition to convection arching and geometry are also important as pointed out by many 
authors in the past. How the large particles move in the exponentially weak convection field 
before being lifted upwards and the 7?,-dependence of the mobility are probably best described 



by Duran et al.'s local arching mechanismsP, 16 . The fact that the whale effect of larger 



particles is less pronounced is certainly due to their lower mobility because of steric hindrance 
effects as formulated by Rosato et al.[^. 

Many of the details of segregation are still not completely clear and in particular in three 
dimensions additional geometrical effects might play a role. This as well as other questions are 
difficult to study conclusively with our numerical technique due to the excessive requirements 
in computer time. It would for instance be interesting to see what happens when the box is so 
wide that the walls of the box are much farther away from the large particle than the height 
of the packing. In this experimentally relevant case the walls would not be able to stabilize 
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Fig. 3. - Strength of the convection rolls measured through the flux $ as a function of the height 
_D, / = 2.6 sec~^ (upper figure), / = 2.8 sec~^ (central figure) and / = 3 sec~^ (lower figure). The 
height D is measured such that the bottom is at the origin of the axis. The flux $ is measured in 
units of particles per period. In the central case the convection rolls are stronger and reach deeper 
inside the material in the presence of the big particle. This triggered convection roll catches the big 
particle and forces it to rise to the top. 



the convection rolls. Simulations with periodic boundary conditions have, however, provided 
rather similar results as with fixed boundaries . It would also be interesting to study larger 
ratios TZ in order to verify predictions made about a characteristic value of 7?. = 12 ||^, ^ 
but for that case one would need to consider substantially larger systems. The limitations in 
observation time due to the computational requirements also puts limits on the determination 
of the segregation velocity and we cannot exclude that particles rise on time scales much larger 
than the ones accessible numerically. 
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Fig. 4. - Evolution of the vertical position of the big particle as function of time, (a) TZ = 2.0 and (b) 
7^ = 3.0 for / = 3.2 sec'^. 
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